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Abstract 

We study analytically and numerically the stability of the standing waves for a non- 
linear Schrodinger equation with a point defect and a power type nonlinearity. A 
main difficulty is to compute the number of negative eigenvalues of the linearized 
operator around the standing waves, and it is overcome by a perturbation method 
and continuation arguments. Among others, in the case of a repulsive defect, we 
show that the standing wave solution is stable in H^^^CK) and unstable in II^{M.) 
under subcritical nonlinearity. Further we investigate the nature of instability: under 
critical or supercritical nonlinear interaction, we prove the instability by blowup in 
the repulsive case by showing a virial theorem and using a minimization method in- 
volving two constraints. In the subcritical radial case, unstable bound states cannot 
collapse, but rather narrow down until they reach the stable regime (a finite-width 
instability). In the non-radial repulsive case, all bound states are unstable, and the 
instability is manifested by a lateral drift away from the defect, sometimes in com- 
bination with a finite-width instability or a blowup instability. 
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1 Introduction 

Solitary waves are localized waves that propagate in nonlinear media where 
dispersion and/or diffraction are present. They appear in various fields of 
physics such as nonlinear optics, Bose- Einstein Condensates (BEC), plasma 
physics, solid state physics, water waves etc. The dynamics of solitons is mod- 
eled by the Nonlinear Schrodinger equation (NLS) in the context of nonlinear 
optics or the Gross-Pitaevskii (GP) equation in the context of BEC. 

By now, the stability and dynamics of solitons in homogeneous media are well 
understood. However, stability and dynamics of solitons in inhomogeneous 
media are still a matter of intense research, both theoretically and experimen- 
tally. Of particular interest is the NLS equation with a linear potential (or 
lattice) 

idtu{t,x) = —d^u — V{x)u — 

(1) 

m(0, x) = Uq. 

In nonlinear optics, the potential V{x) describes the variation of the linear 
refractive index in space. In BEC, it describes an external potential applied 
to the condensate. The potential can be localized (e.g., a single waveguide 
in nonlinear optics [39,32]), parabohc (e.g., a magnetic trap in BEC [1,31]) 
or periodic (e.g., a waveguide array or photonic crystal lattice in nonlinear 
optics [42]). 

In the presence of a potential, a key parameter is the relative width of the 
solitary wave, compared with the characteristic length-scale of the potential. 
For example, in the case of a periodic lattice, narrow solitary waves are af- 
fected, to leading order, by the local changes of the potential near the soliton 
center [16,14,15,43], whereas wide solitary waves are affected by the potential 
average over a single period [14,15]. 

In this paper we consider a NLS/GP equation with a delta function potential 

idtu{t, x) = —diu — '-fS{x)u — \u\P~^u, 

(2) 

u{0, x) = Uq, 
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where 7 G M, 1 < p < +00 and (t, x) G M"*" x M. Here, S is the Dirac 
distribution at the origin, namely, {S,v) = v{0) for v G if^(]R). Equation (2) 
can be viewed as a prototype model for the interaction of a wide soliton with 
a narrow potential. The main advantage of using the delta-function potential 
rather than a finite-width potential is the existence of an explicit expression for 
the soliton profile. This allows to prove results, whose proofs are considerably 
harder for a general linear potential. 

In nonlinear optics, Eq. (2) models a soliton propagating in a medium with a 
point defect [22,36] or the interaction of a wide soliton with a much narrower 
one in a bimodal fiber [6]. In BEC, this equation models the dynamics of a 
condensate in the presence of an impurity of a length-scale much smaller than 
the healing length. Such an impurity can be realized by a tightly focused beam, 
by another spin state of the same atom or by another alkali atom confined in an 
optical trap [41]. In contrast to wide solitons in a periodic potential, in Eq. (2) 
the (wide) soliton profile is affected only by the local variation of the potential 
rather than by its average. Moreover, since the potential is localized, there is 
no band structure and gap solitons characteristic of a periodic potential, see 
e.g., [10]. 

Equation (2) was studied previously by several authors. 
In [6,22,25,26,27,40,41], the phenomenon of soliton scattering by the ef- 
fect of the defect was observed, namely, interactions between the defect 
and the homogeneous medium soliton. For example, varying amplitude and 
velocity of the soliton, they studied how the defect is separating the soliton 
into two parts : one part is transmitted past the defect, the other one is 
captured at the defect. Holmer, Marzuola and Zworski [25,26] gave numerical 
simulations and theoretical arguments on this subject. Recently, these results 
were observed experimentally for a single waveguide potential [32]. 

In this paper, we study the stability and instability of the standing-wave so- 
lution of (2) of the form u{t,x) = e^^^p{x) where ip is required to satisfy 

^Gi/i(K)\{0}. 

V 

Stability under radial (symmetric) perturbations was studied analytically 
in [22,13,12]. In this paper, we study stability under non-ra(iza/ perturbations. 
We also show that the instability associated with momentum-nonconserving 
perturbations is excited only for a repulsive defect (7 < 0), and is manifested 
by a lateral movement of the wave away from the defect. 

In the numerical part of this study we combine some recent ideas such as a 
quantitative approach to (in) stability and characterization of the instability 
type (width or drift instability) in order to provide a systematic description 
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of the standing wave dynamics. We emphasize that both our approach and 
results are relevant to standing waves of the NLS (1) with a general linear 
potential, and also to NLS with a nonlinear potential [14,15]. 



2 Review of previous results 



Notations: The space L''(]R, C) will be denoted by L^(R) and its norm by 
II ■ \\r. When r = 2, the space L^(]R) will be endowed with the scalar product 



{u,v)2 = Re / uvdx for u,v G L^(M). 

JR 



The space H^(M., C) will be denoted by H^{R), its norm by || ■ ||_h-i(r) and the 
duality product between H-^{R) and H^{R) by (-,■). We write H^^^{R) for 
the space of radial (even) functions of (M) : 

H^^^{R) = {v e H\R); v{x) = v{-x), x G R}. 



When 7 = 0, the set of solutions of (3) has been known for a long time. 
In particular, modulo translation and phase, there exists a unique positive 
solution, which is explicitly known. This solution is even and is a ground state 
(see, for example, [4,7,29] for such results). When 7 7^ 0, an explicit solution 
of (3) was presented in [13,22] and the following was proved in [12,13]. 

Proposition 1 Let uj > 7^/4. Then there exists a unique positive solution 
V'a;,7 of (3) . This solution is the unique positive minimizer of 



■ml{S^^,{v)- V e H\R) \ {0}, I^^,{v) = 0} 7 ^ 0, 

inf {S^,,{v)- V E \ {0}, L,,{v) = 0} 7 < 0, 



d{uj) = 

where S^^-y and 1^^^ are defined for v E (M) by 



1 no ..112 , ^11 ||2 7, /nM2 1 II iip+l 



Suniv) = i^\\dxv\\2 + -\\v\\2 - 2 K'(o)l - ^:^ll^'lip+i' 

luni'") = \\dxv\\l + t^||f II2 - 71^(0)1^ - IK'llp+i- 



Furthermore, we have an explicit formula for if, 



i^±l)^sech^fi^^^ll^H+tanh-^ ^ 



2^^7 



1 

p-1 



(4) 
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The dependence of ^9^^^ on uj and 7 can be seen in Figure 1. The parameter uj 
affects the width and height of v?(^,7: the larger lo is, the narrower and higher 
Lp^^^^ becomes, and vise versa. The sign of 7 determines the profile of (/?(^^ near 
x = 0: It has a "V" shape when 7 < 0, and a "A" shape when 7 > 0. 



(b) 
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Fie. 1. function of x for a; = 4 (solid line) and uj = 0.5 (dashed line), (a) 

7 = 1; (b) 7 = —1. Here, p = 4. 

Remark 1 (i) As it was stated in [12, Remark 8 and Lemma 26], the set 
of solutions of (3) 

{v e H^{M) \ {0} such that - dlv + uv - -^v5 - \v\^~^v = 0} 

is explicitly given by {e^^ipu},-/] ^ G M}. 
(ii) There is no nontrivial solution in if^(M) for u ^ 7^/4. 

The local well-posedness of the Cauchy problem for (2) is ensured by [7, The- 
orem 4.6.1]. Indeed, the operator —d"^ — 7^ is a self-adjoint operator on L^(M) 
(see [2, Chapter 1.3.1] and Section 2 for details). Precisely, we have 

Proposition 2 For any Uq G H^(R), there exist T^g > and a unique 
solution u e C{[0,T^^),H\R)) n C\[0,T^^), H-\R)) of (2) such that 
limf^TuQ \\dxu\\2 = +C)0 if < +00. Furthermore, the conservation of en- 
ergy and charge hold, that is, for any t G [0,T„g) we have 



E{u{t)) = E{uo), (5) 
Ht)\\l = ho\\l, (6) 

where the energy E is defined by 

i?(^) = ^l|5.^ll^-|l^(0)p-^||^^|lSl, for veH\R). 



(see also a verification of this proposition in [13, Proposition 1]). 

Remark 2 From the uniqueness result of Proposition 2 it follows that if an 
initial data Uq belongs to H^^^(M.) then u(t) also belongs to H^^^{M.) for all 
tG [0,T„J. 
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We consider the stability in the following sense. 

Definition 3 Let ip be a solution of (3). We say that the standing wave 
u{x,t) = e'"^'^ip{x) is (orbitally) stable in H^(M.) (resp. H^^^(M.)) if for any 
e > there exists r] > with the following property : if Uo G if^(M) (resp. 
H^^^{R)) satisfies \\uo — </?||_f/i(]R) < rj, then the solution u{t) of (2) with 
u(0) = Uq exists for any t > and 

sup inf \\u{t) — e'^(/?||iyi(]R) < e. 
te[o,+oo) 

Otherwise, the standing wave u{x,t) = e*'^V(^) ^■^ said to be (orbitally) un- 
stable m H\R) (resp. H^^CR)). 

Remark 4 With this definition and Remark 2, it is clear that stability in 
H^{R) implies stability in H^g_^{R) and conversely that instability in H^.^^{R) 
implies instability in H^(R). 

When 7 = 0, the orbital stability for (2) has been extensively studied (see 
[3,7,8,44,45] and the references therein). In particular, from [8] we know that 
e*'^Vaj,o(3^) is stable in if^(R) for any Ci;>0ifl<|)<5. On the other hand, 
it was shown that e*'^Va;,o(2^) is unstable in H^{M.) for any > if j> ^ 5 (see 
[3] for p > 5 and [45] for p = 5). 

In [22], Goodman, Holmes and Weinstein focused on the special case p = 3, 
7 > and proved that the standing wave e*'^*(y9t<j,7(^) is orbitally stable in 
H^{R). When 7 > 0, the orbital stability and instability were completely 
studied in [13] : the standing wave e*'^V(.^,7(^) is stable in if^(M) for any 
(J > 7^/4 if 1 < p ^ 5, and if p > 5, there exists a critical frequency uji > 7^/4 
such that e*'^Vw,7(3^) is stable in H^(W) for any u G (7^/4, a;i) and unstable 
in H^{M.) for any uj > Ui. 

When 7 < 0, Fukuizumi and Jeanjean showed the following result in [12]. 

Proposition 3 Let 7 < and uj > 7^/4. 

(i) If 1 < p ^ 3 the standing wave e*'^V<.^,7(^) ^■^ stable in H^^^(M.). 
(a) If 3 < p < 5, there exists 002 > 7^/4 such that the standing wave 
G^^^^Lorii^) "^s stable in H^.^^{W) when uj > UJ2 and unstable in iJ^(]R) when 
7^/4 <u} <u]2- 

(Hi) If p ^ 5, then the standing wave e^^^ip^^^{x) is unstable in i/^(]R). 
The critical frequency UJ2 is given by 
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J(w2)(p-5) 7 



P 



71 



-(p-3)/(p-l) 



+00 



sec 



h^'^''~^\y)dy, A{uj2,i) = tank"' 
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3 Summary of results 



The results of stability of [12] recalled in Proposition 3 assert only on stabil- 
ity under radial perturbations. Furthermore, the nature of instability is not 
revealed. In this paper, we prove that there is instability in the whole space 
when stability holds under radial perturbation (see Theorem 4), and that, 
when p ^ 5, the instability established in [12] is strong instability (see Defini- 
tion 6 and Theorem 5). 

Our first main result is the following. 

Theorem 4 Let 7 < and u > 7^/4. 

(i) If 1 < p ^ 3 the standing wave e*'^*v?i^,7(^) ^■^ unstable in iJ^(]R). 
(a) If 3 < p < 5, the standing wave e*'^V<^,7(^) ^■^ unstable in H^(R) for any 
uj > UJ2, where 002 is defined in Proposition 3. 

As in [12,13], our stability analysis relies on the abstract theory by Grillakis, 
Shatah and Strauss [23,24] for a Hamiltonian system which is invariant under 
a one-parameter group of operators. In trying to follow this approach the main 
point is to check the following two conditions: 

(1) The slope condition: The sign of 9a;||v9^,^||2. 

(2) The spectral condition: The number of negative eigenvalues of the lin- 
earized operator 

Llu;V = -dlv + UJV- -idv - p^pI^^v. 

We refer the reader to Section 4 for the precise criterion and a detailed ex- 
planation on how Lj^ appears in this stability analysis. Making use of the 
exphcit form (4) for (pu},-y, the sign of 9(^||v5w,7ll2 "was explicitly computed in 
[12,13]. 

In [12], a spectral analysis is performed to count the number of negative eigen- 
values, and it is proved that the number of negative eigenvalues of Lj^^ in 
H^^^iM.) is one. This spectral analysis of Lj^^^ is relying on the variational 
characterization of (fuj.-y- However, since ^9^^^ is a minimizer only in the space 
of radial (even) functions H^^{R), the result on the spectrum holds only in 
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H}^^(M.), namely for even eigenf unctions. Therefore the number of negatives 
eigenvalues is known only for Lj^ considered in H^^^{M.). With this approach, 
it is not possible to see whether other negative eigenvalues appear when the 
problem is considered on the whole space H^{R). 

To overcome this difficulty, we develop a perturbation method. In the case 
7 = 0, the spectrum of L^^ is well known by the work of Weinstein [46] (see 
Lemma 14) : there is only one negative eigenvalue, and is a simple isolated 
eigenvalue (to see that, one proves that the kernel of L^^ is spanned by d^^Puifiy 
that i9x</5a;,o has only one zero, and apply the Sturm Oscillation Theorem). 
When 7 is small, Lj^^ can be considered as a holomorphic perturbation of 
L5^. Using the theory of holomorphic perturbations for linear operators, we 
prove that the spectrum of Lj^ depends holomorphically on the spectrum 
of 

^luj Lemma 15). Then the use of Taylor expansion for the second 
eigenvalue of Lj^^ allows us to get the sign of the second eigenvalue when 7 
is small (see Lemma 16). A continuity argument combined with the fact that 
if 7 7^ the nullspace of Lj^ is zero extends the result to all 7 G M (see the 
proof of Lemma 12). See subsection 4.2 for details. We will see that there are 
two negative eigenvalues of Lj^ in if ^(M) if 7 < 0. 

Remark 5 (i) Our method can be applied as well in ii^(M) or in H}^^{W), 
and for 7 negative or positive (see subsections 4.3 and 4.4). Thus we 
can give another proof of the result of [13] in the case 7 > and of 
Proposition 3. 

(ii) The study of the spectrum of linearized operators is often a central point 
when one wants to use the abstract theory of [23,24]. See [14,17,18,19,28] 
among many others for related results. 

The results of instability given in Theorem 4 and Proposition 3 say only that 
a certain solution which starts close to (fcu,-y will exit from a tubular neighbor- 
hood of the orbit of the standing wave in finite time. However, as this might be 
of importance for the applications, we want to understand further the nature 
of instability. For that, we recall the concept of strong instability. 

Definition 6 A standing wave e^'^^ip{x) of (2) is said to he strongly unstable 
in H^{M) if for any e > there exist G H^{M.) with \\ue — ^\\h'^{K) < ^ o-^^d 
< +00 such that Xim^Tu^ \\dxu{t)\\2 = +00, where u(t) is the solution of 
(2) with u{0) = Ue- 

Our second main result is the following. 

Theorem 5 Let 7 ^ 0, > 7^/4 and p ^ 5. Then the standing wave 
e*'^V(^.7(^) ^■^ strongly unstable in H^(R). 

Whether the perturbed standing wave blows-up or not depends on the pertur- 
bation. Indeed, in Remark 10 we define an invariant set of solutions and show 



8 



that if we consider an initial data in this set, then the solution exists globally 
even when the standing wave e*'^*</)(^ is strongly unstable. 

We also point out that when 1 < p < 5, it is easy to prove using the con- 
servation laws and Gagliardo-Nirenberg inequality that the Cauchy problem 
in if^(]R) associated with (2) is globally well posed. Accordingly, even if the 
standing wave may be unstable when 1 < p < 5 (see Theorem 4), a strong 
instability cannot occur. 

As in [3,45], which deal with the classical case 7 = 0, we use the virial identity 
for the proof of Theorem 5. However, even if the formal calculations are similar 
to those of the case 7 = 0, a rigorous proof of the virial theorem does not 
immediately follow from the approximation by regular solutions (e.g. see [7, 
Proposition 6.4.2], or [20]). Indeed, the argument in [7] relies on the if^(M) 
regularity of the solutions of (2). Because of the defect term, we do not know if 
this if^(M) regularity still hold when 7 7^ 0. Thus we need another approach. 
We approximate the solutions of (2) by solutions of the same equation where 
the defect is approximated by a Gaussian potential for which it is easy to have 
the virial theorem. Then we pass to the limit in the virial identity to obtain : 

Proposition 6 Let uq G H^{M.) such that xuq G L^(M) and u(t) be the solu- 
tion of (2). Then the function f : t \\xu(t)\\2 is C"^ and 



dtf (t) = Aim / uxdxudx, (7) 
d^f{t)=8Q,{u{t)), (8) 

where is defined for v G H^{M.) by 



Even if we benefit from the virial identity, the proofs given in [3,45] for the case 
7 = do not apply to the case 7 < 0. For example, the method of Weinstein 
[45] in the case p = 5 requires in a crucial way an equality between 2E and Q 
which does not hold anymore when 7 < 0. Moreover, the heart of the proof 
of [3] consists in minimizing the functional Si^^^ on the constraint Q^yiv) = 0, 
but the standard variational methods to prove such results are not so easily 
applied to the case of 7 7^ 0. To get over these difficulties we introduce an ap- 
proach based on a minimization problem involving two constraints. Using this 
minimization problem, we identify some invariant properties under the flow of 
(2). The combination with these invariant properties and the conservation of 
energy and charge allows us to prove strong instability. We mention that some 
related techniques have been introduced in [33,34,35,37,47]. In conclusion, we 
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can give a simpler method to prove Theorem 5 than that of [3] even though 
we have a term of delta potential. 

Remark 7 The case 7 < 0, = u;2 and 3 < p < 5 cannot be treated with our 
approach and is left open (see Remark 8). In light of Theorem 4, we believe 
that the standing wave is unstable in this case, at least in H^{M.) (see also [12, 
Remark 12]). When 7 > 0, the case lj = ui and p > 5 is also open (see [13, 
Remark 1.5]). 

Let us summarize the previously known and our new rigorous results on sta- 
bility in (2): 

(i) For both positive and negative 7, there is always only one negative eigen- 
value of linearized operator in H^^^{R) ([12], subsection 2.5). Hence, the 
standing wave is stable in H^^^{R) if the slope is positive, and unstable 
if the slope is negative. 

(ii) 7 > 0. In this case the number of the negative eigenvalues of linearized 
operator is always one in H^{M.). Stability is determined by the slope 
condition, and the standing wave is stable in H^^^(M.) if and only if it is 
stable in H^iR). Specifically ([12,13], subsection 2.4), 

(a) 1 < p < 5: Stability in H^{R) for any u > 7^/4. 

(b) 5 <p: Stability in H^{R) for 7^/4 < cu < toi, instability in H^^^(R) 
for uj > ui. 

(iii) 7 < 0. In this case the number of negative eigenvalues is always two 
(Lemma 12) and all standing waves are unstable in H^{R) (Theorem 4 
and Theorem 5). Stability in if^\j(M) is determined by the slope condi- 
tion and is as follows ([12]): 

(a) 1 < p < 3: Stability in H^^^(R) for any u > 774. 

(b) 3 < p < 5: Stability in H^^^(M.) for u > UJ2, instability in H^^^{R) for 
7^/4 < UJ < UJ2. 

(c) 5 < p: Strong instability in H^^^{R) (and in H^{M.)) for any 7^/4 < 
(Theorem 5). 

There are, however, several important questions which are still open, and 
which we explore using numerical simulations. Our simulations suggest the 
following: 

(i) Although an attractive defect (7 > 0) stabilizes the standing waves in 
the critical case {p = 5), their stability is weaker than in the subcritical 
case, in particular for < 7 <^ 1. 

(ii) Theorem 5 shows that instability occurs by blowup when 7 < and 
p ^ 5. In all other cases, however, it remains to understand the nature 
of instability. Our simulations suggest the following: 

(a) When 7 > 0, p > 5, and u > ui, instability can occur by blowup. 

(b) When 7 < 0, 3 < p < 5, and 7^/4 < u < U2, the instability in 



10 



i74,d(^) is a. finite-width instability, i.e., the solution initially narrows 
down along a curve 4'uj*{t),'y, where uj*{t) can be defined by the relation 

max(/)^.(t),^(x) = max\u{x,t)\. 

As the solution narrows down, w*{t) increases and crosses from the 
unstable region u; < ti;2 to the stable region ou > lj2- Subsequently, 
collapse is arrested at some finite width, 
(c) When 7 < 0, the standing waves undergo a drift instability, away 
from the (repulsive) defect, sometimes in combination with finite- 
width or blowup instability. Specifically, 

(c.i) When 1 < p < 3 and when 3 < p < 5 and u > UJ2 (i.e., when 
the standing waves are stable in H^^^(M.)), the standing waves 
undergo a drift instability. 
(c.ii) When 3 < p < 5 and 7^/4 < u < UJ2, the instability in H^{R) is 
a combination of a drift instability and a finite- width instability, 
(c.iii) When p > 5, the instability in H^{M.) is a combination of a drift 
instability and a blowup instability, 
(iii) Although when p = 5 and 7 > 0, and when p > 5, 7 > 0, and 
7^/4 < u < uji the standing wave is stable, it can collapse under a 
sufficiently large perturbation. 

We note that all of the above holds, more generally, for NLS equations with a 
nonlinear potential [14,15] and for narrow solitons of a linear potential [43]. 

The rest of the paper is organized as follows. In Section 4, we prove Theorem 
4 and explain how our method allows us to recover the results of [12,13]. In 
Section 5, we establish Theorem 5. Numerical results are given in Section 6. 

Throughout the paper the letter C will denote various positive constants whose 
exact values may change from line to line but are not essential to the analysis 
of the problem. 
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4 Instability with respect to non-radial perturbations 

We use the general theory of Grillakis, Shatah and Strauss [24] to prove The- 
orem 4. 

First, we explain how we derive a criterion for stability or instability for our 
case from the theory of Grillakis, Shatah and Strauss. In our case, it is clear 
that Assumption 1 and Assumption 2 of [24] are satisfied. The last assump- 
tion. Assumption 3, should be checked. We consider the sesquilinear form 
S'^,'yi^uj,y) ■■ H\R) X H\R) ^ C as a linear operator : H\R) H-\R). 
The spectrum of is the set {A G C such that H[l — XI is not invertible}, 
where / denote the usual if^(]R) — if~^(M) isomorphism, and we denote by 
n{H2) the number of negative eigenvalues of H^- Having established the as- 
sumptions of [24], the next proposition follows from [24, Instability Theorem 
and Stability Theorem]. 

Proposition 7 (1) The standing wave e*'^"*v9(^o,7(^) ^■^ unstable if the integer 
{n{H2 ) —p{d"{uJo)) is odd, where 



(2) The standing wave e*'^°VaJo,7(^) ^■^ stable if {n{H2^) — p{d" {ujq)) = 0. 

Let us now consider the case 7 < 0. It was proved in [12] that 

Lemma 8 Let 7 < and uj > 7^/4. We have : 

(i) // 1 < p ^ 3 and u > 7^/4 then d^\\ip^^^\\l > 0, 

(a) If 3 < p < 5 and u > 002 then d^\\ip^^^\\2 > 0, 

{Hi) If 3 < p < 5 and 7^/4 < uj < LJ2 then duj\\(puj,'y\\2 < 0, 

(iv) If p ^ 5 and u > 7^/4 then d:^\\ipuj,-y\\l < 0. 

Thus Theorem 4 follows from Proposition 7, Lemma 8 and 

Lemma 9 1/7 < 0, then n{H2) = 2. 

Remark 8 (1) Let 7 < 0. In the cases 3 < p < 5 and u; < ci;2 or p > 5 it 
was proved in [12] that i9tj||9?aj,7ll2 < 0- From Lemma 9, we know that the 
number of negative eigenvalues of is n{H2) = 2 when is considered 
on the whole space if^(]R). Therefore n{H^) — p{d"{uj)) = 2 and this 
corresponds to where the assumption of [24] may not be applied. 

However, if we consider in II^^^{M.), then it follows from [12] that 
n{H2) = 1, thus n{H2) —p{d"{uj)) = 1. Then, we can apply Proposition 
7 to this case and it allows us to conclude instability in II^^^{M.) (as it 
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was done in [12]). But, with Remark 4, we can conclude that instability 
holds on the whole space if^(]R). 
(2) Note that the case uj = ijJ2 corresponds to 9i^||v5(^,7ll2 = (3 < p < 5) 
and will not be treated here. In view of Theorem 4, we believe that the 
standing wave is unstable in this case, at least in iJ^(M). 

We divide the rest of this section into four parts. In subsection 4.1 we introduce 
the general setting to perform our proof, and we study whether Assumption 
3 of [24] is satisfied. Lemma 9 will be proved in subsection 4.2. Finally, we 
discuss the positive case and the radial case in subsections 4.3 and 4.4. 

4-1 Setting for the spectral problem 

To express H^, it is convenient to split u in real and imaginary part : for 
u G H^{R,C) we write u = ui + iu2 where ui = Re{u) e if-'^(M, M) and 
U2 = lm{u) E H\R,R). Now we set 

where the operators L\^, : H^{'K, M) ^ H'^{W) are defined for v e H^{R) 
by 

Ll^^v = -dlv + UJV- 'yvS - p^l'^v, 

When we will work with Lj^^, L^^^, the functions considered will be understood 
to be real valued. 

For the spectral study of H^, it is convenient to view as an unbounded 
operator on L^(M), thus we rewrite our spectral problem in this setting. First, 
we redefine the two operators Lj^ and as unbounded operators on L^(M). 
We begin by considering the bilinear forms on H^(M.) associated with Lj^ and 
-^2,w by setting for v, w G H^{M.) 

Bi^u^{v,w) := (l1^^v,w) and Bl^{v,w) := iL\^v,w) , 
which are explicitly given by 

Bi^ujiv.w) = J^d^vd^wdx+ujJ^vwdx--fv{Q)w{Q)-J^p(pP-^vwdx, 

[9) 

Bloj{v,w) = jji^d^vd:,wdx+ujjji^vwdx--fv{0)w{0) - Jji^(pP~^^vwdx. 

Let us now consider Bj^ and i^Jcj as bilinear forms on L^(M) with domain 
D{Bl^^) = DiyBl^^) := if^(M). It is clear that theses forms are bounded 
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from below and closed. Then the theory of representation of forms by op- 
erators (see [30, VI. §2.1]) implies that we define two self-adjoint operators 
■■ Di^^J C L\R) ^ L\R) and lJ, : D(lJ'J C L\R) L\R) by 
setting 

D{Ll^):={v e H\R)\3w G L\R) s.t. Wz G H\R), Bl^{v,z) = {w,z)2}, 
Dil^J ■■={v G H\R)\3w G L\R) s.t. \fz G H\R), Bl^{v,z) = {w.z)^}. 

and setting for v G Di^Lj^^) (resp. v G D^Ll^^^)) that L\^^v := (resp. 

-^2,w^ •= where w is the (unique) function of L^(]R) which satisfies 
Bl^iv.z) = {w,z)2 (resp. Bl^{v,z) = {w,z)2) for all z G H^R). 

For notational simplicity, we drop the tilde over Lj^ and L2^. It turns out 
that we are able to describe explicitly Lj^ and -^2^. 

Lemma 10 The domain of Lj^^ and of L^i^ in L^(M) is 

D^ = {ve H\R) n H'^iR \ {0}); 9^t;(0+) - (9^t;(0") = -7t;(0)} 
anc? /or f G the operators are given by 

LlujV = -dlv + uv - p^l'^v, ^^^^ 
Llu^v = -dlv + ujv - '^P-^v. 

The proof of this lemma is given in Appendix A. We conclude this subsection 
mentioning some basic properties of the spectrum of H2- Precisely, to check 
[24, Assumption 3] is equivalent to check the following lemma. 

Lemma 11 Lei 7 G M \ {0} and uj > 7^/4. 

(i) The operator H^, has only a finite number of negative eigenvalues, 

(ii) The kernel of is span{i(p^ .y} , 

(iii) The rest of the spectrum of is positive and bounded away from 0. 

Our proof of Lemma 11 borrows some elements of [12]. In particular, (ii) 
in Lemma 11 is shown in [12, Lemma 28 and Lemma 31]. For the sake of 
completeness, we provide a proof in Appendix B. 

4-2 Count of the number of negative eigenvalues 

In this subsection, we prove Lemma 9. First, we remark that, as it was shown in 
the proof of Lemma 11, is the first eigenvalue of I/2,a;. Thus n{H2) = n^Lj^^), 
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where n{Ll^J) is the number of negative eigenvalues of L\^^. Therefore, Lemma 
9 follows from 

Lemma 12 Let 7 < and u > 7^/4. Then n{L'l^^) = 2. 

Our proof of Lemma 12 is divided in two steps. First, we use a perturbative 
approach to prove that, if 7 is close to and negative, L^^ has two negative 
eigenvalues (Lemma 16). To do this, we have to ensure that the eigenvalues and 
the eigenvectors are regular enough with respect to 7 (Lemma 15) to make use 
of Taylor formula. It follows from the use of the analytic perturbation theory 
of operators (see [30,38]). The second step consists in extending the result of 
the first step to any values of 7 < 0. Our argument relies on the continuity of 
the spectral projections with respect to 7 and it is crucial, as it was proved 
in Lemma 11, that can not be an eigenvalue of Lj^ (see [17,18] for related 
arguments) . 

We fix > 7^/4. For the sake of simplicity we denote Lj^^ by Lj and (puj^-y by 
ip^, and so on in this section 2. The following lemma verifies the holomorphicity 
of the operator Lj^, see proof in Appendix B. 

Lemma 13 ^45 a function of 7, (Lj) is a real-holomorphic family of self- 
adjoint operators (of type (B) in the sense of Kato). 

The following classical result of Weinstein [46] gives a precise description of 
the spectrum of the operator we want to perturb. 

Lemma 14 The operator L\ has exactly one negative simple isolated first 
eigenvalue. The second eigenvalue is 0, and it is simple and isolated. The 
nullspace is span{c}a;(y9o} , and the rest of the spectrum is positive. 

Combining Lemma 13 and Lemma 14, we can apply the theory of analytic per- 
turbations for linear operators (see [30, VII. §1.3]) to get the following lemma. 
Actually, the perturbed eigenvalues are holomorphic since they are simple. 

Lemma 15 There exist 70 > and two functions A : (—70, 7o) ^ (^i^d 
f : (—70,70) ^ L^(M) such that 

(i) A(0) = and /(O) = d^ipo, 

(ii) For all 7 G (— 7o,7o)j -^(7) is the simple isolated second eigenvalue of Lj 
and f{'-f) is an associated eigenvector, 

(iii) A(7) and f{'-f) are holomorphic in (— 7o,7o)- 

Furthermore, 70 > can be chosen small enough to ensure that, expect the 
two first eigenvalues, the spectrum of Lj is positive. 
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Now we investigate how the perturbed second eigenvalue moves depending on 
the sign of 7. 

Lemma 16 There exists < 71 < 70 such that A(7) < for any —71 < 7 < 
and A(7) > for any < 7 < 71. 

Proof of Lemma 16. We develop the functions A(7) and /(7) of Lemma 15. 
There exist Aq G M and /o G L^(]R) such that for 7 close to we have 

A(7)=7Ao + 0(7'), (11) 
/(7) = 9x^0 + 7/0 + 0(7'). (12) 

From the explicit expression (4) of y?^, we deduce that there exists G H^{M) 
such that for 7 close to we have 

(yj^ = (y9o + 75-0 + 0(7^). (13) 

Furthermore, using (13) to substitute into (3) and differentiating (3) with 
respect to 7, we obtain 

(l?^7o,V') = V'o(0)^(0) + 0(7), (14) 

for any ^jj G iJ^(M). 

To develop Aq with respect to 7, we compute {Lif{'y),dx(po)2 in two different 
ways. 

On one hand, using Ljf^'j) = A(7)/(7), (11) and (12) lead us to 

(Ll/(7), d,ipo)2 = XolWd.Ml + 0(f). (15) 

On the other hand, since Lj is self-adjoint, we get 

mf{^),d,ipo)2 = {f{l),Ljd,^o)2. (16) 

Here we note that d,j:(po G D{L\) : indeed, d^^Q G H'^iM) and dx^Q{{i) = 0. 
We compute the right hand side of (16). We use (10), L^d^^Po = 0, and (13) 
to obtain 

= -jp{p - l)^l'god,^o + 0(7'). (17) 
Hence, it follows from (12) that 

mfil), dx^o)2 = -{dx<Po, igop{p - l)^l~'^dx^o)2 + 0(7^). (18) 
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Now, as it was remarked in [14, Lemma 28], it is easy to see that using (3) 
with 7 = we get 



LiiVo -Vo ^) = P{P - l)v'o ^^xV'o' 



(19) 



which combined with (18) gives 



(L7/(7),a,.^o)2 = -7(^' 



'^o,<^o-<^o) + C'(7^). 



(20) 



Finally, with (14) we obtain from (20) 



7(v.o(0)^- ^0(0)^+^) + 0(f). 



(21) 



Combining (21) and (15) we obtain 



,p+i 



+ 0(7). 



ll'^xV'olli 



It follows that Aq is positive for sufficiently small I7I, which in view of (11) 
ends the proof. □ 

We are now in position to prove Lemma 12. 
Proof of Lemma 12. Let 700 be defined by 

7oo = inf{7 < 0; has exactly two negative eigenvalues for all 7 G (7,0]}. 

From Lemma 16, we know that 700 is well defined and 700 G [—00, 0). Arguing 
by contradiction, we suppose 700 > —00. 

Let N be the number of negative eigenvalues of 1^1°° . Denote the first eigen- 
value of ^1°° by A^^. Let F be defined by 

F = {2; G C; z = Zi+iz2, {zi, Z2) G [—6, 0]x [—a, a], for some a > 0, 6 > |A^_^|}. 

From Lemma 11, we know that L\°^ does not admit zero as eigenvalue. Thus F 
define a contour in C of the segment [A^^ , 0] containing no positive part of the 
spectrum of ^1°° , and without any intersection with the spectrum of Vl°° . It is 
easily seen (for example, along the lines of the proof of [30, Theorem VII-1.7]) 
that there exists a small 7* > such that for any 7 G [700 — 7*, 7oo + 7*], we 
can define a holomorphic projection on the negative part of the spectrum of 
Li contained in F by 



Let us insist on the fact that we can choose F independently of the parameter 
7 because is not an eigenvalue of LI for all 7. 
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Since 11 is holomorphic, 11 is continuous in 7, then by a classical con- 
nectedness argument (for example, see [30, Lemma 1-4.10]), we know that 
dim(Ran 11(7)) = N for any 7 G [700 — 7*,7oo + 7*]- Furthermore, is ex- 
actly the number of negative eigenvalues of Lj when 7 G [700 — 7*,7oo + 7*] 
: indeed, if Lj has a negative eigenvalue outside of F it suffice to enlarge F 
(i.e., enlarge b) until it contains this eigenvalue to raise a contradiction since 
then L7°° would have, at least, iV-fl eigenvalues. Now by the definition of 700, 
^7oo+7* ]-^g^g negative eigenvalues and thus we see that Lj has two negative 
eigenvalues for all 7 G [700 — 7*, 0[ contradicting the definition of 700- 

Therefore 700 = — 00. □ 

Remark 9 In [12, Lemma 32], the authors proved that there are at most two 
negative eigenvalues of Lj in H^{M.) using variational methods. In our present 
proof, we can directly show that there are exactly two negative eigenvalues 
without such variational techniques. 

4-3 The case 7 > 

The proof of Lemma 12 can be easily adapted to the case 7 > 0, and with 
Lemma 16 we can infer that Lj has only one simple negative eigenvalue when 
7 > 0. Since n{H"') = n(L]'), it follows that (in the following Lemmas 17, 18 
and Proposition 19, there is no omission of parameter u to understand the 
dependence clearly) 

Lemma 17 Let 7 > and u > 7^/4. Then the operator has only one 
negative eigenvalue, that is n{H^) = 1. 

When 7 > 0, the sign of 9^||v:'a;,7ll2 computed in [13]. Precisely : 
Lemma 18 Let 7 > and uj > 7^/4. We have : 

(i) // 1 < p ^ 5 and u > 7^/4 then d^\\ip^^y\\l > 0, 
(a) If p > 5 and 7^/4 < uj < Ui then d^\\ip^^y\\2 > 0, 
{Hi) If p > 5 and uj > Ui then duj\\(puj,'y\\2 < 0. 

Here uji is defined as follows: 



p — 5 



7 



( 




) 



(p-3)/{p-l) 



P — 1 



J{uJi) 




1 



JM = f 

J A 



00 



A(t^i,7) 
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Then, using Lemma 17, Lemma 18 and Proposition 7, we can give an alterna- 
tive proof of [13, Theorem 1] (see also [12, Remark 33]). Precisely, we obtain: 

Proposition 19 Let 7 > and tu > 7^/4. 

(i) Let 1 < p ^ 5. Then e*'^Vw,7(^) ^■^ stable in H^{M.) for any uj G (7^/4, +00). 
(m) Let p > 5. Then e^'^''^ip^^^{x) is stable in H^(M.) for any u G (7^/4, cji), and 
unstable in if^(]R) for any uj G (cji, +00). 



4.4 The radial case 



Before we start to discuss the stability in the radial case, we mention the 
following remarkable fact. 

Lemma 20 The function 7(7) defined in Lemma 15 and corresponding to the 
second negative eigenvalue of Li can be extended to (—00, +00). Furthermore, 
f{'~f) G H^{M.) is an odd function, for each 7 G (—00, +00). 

The proof uses a similar idea to that of Lemma 12, see Appendix C. 

We can deduce the number of negative eigenvalues of Lj in the radial case 
from the result on the eigenvalues of Lj considered in the whole space L^(M). 
Indeed, Lemma 20 ensures that the second eigenvalue of Lj considered in the 
whole space L^(M) is associated with an odd eigenvector, and thus disappears 
when the problem is restricted to subspace of radial functions. Furthermore, 
since ip^ G H}^^{W) and (Ljip^jip^) < 0, we can infer that the first negative 
eigenvalue of Lj is still present when the problem is restricted to sets of radial 
functions. Recalling that n{H^) = n(L]'), we obtain. 

Lemma 21 Let 7 < 0. Then the operator H'' considered on H^^^{M.) has only 
one negative eigenvalue, that is n{H"') = 1. 

Combining Lemma 21, Lemma 8 and Proposition 7, we recover the results of 
[12] recalled in Proposition 3. 

Alternatively, subsection 4.2 can be adapted to the radial case. All the function 
spaces should be reduced to spaces of even functions, and Lemma 21 can also 
be proved in this way. 
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5 Strong instability 



This section is devoted to the proof of Theorem 5. We use the virial theorem 
(Proposition 6) whose verification will be given in Appendix D. 

We begin by introducing some notations 

^ = {t; G Hl^{R) \ {0}; Q^{v) = 0, L,,{v) ^ 0}, 

where 5*^^^^ and I^^^y are defined in Proposition 1 and in Proposition 6. 
Our proof is divided in three steps. 

Step 1. We prove that (f^^^^ is also a minimizer of d^. 

Because of Pohozaev identity Q'y{fuj,'y) = (see [4]), it is clear that ^ '^(i^), 
thus we only have to show ^ d^u). Let v G If Iuj,-y{v) = 0, we have 
Suj,^{v) ^ d{Lj), therefore we suppose /(^.^(f) < 0. For a > 0, let v"" be such 
that f"(x) = a^/^f (ax). We have 

Lni^n = a'\\d,v\\l + uj\\v\\l - ^a\vm' - "^""'^^I^IlK' 

thus lim Iuj-^f(v") = (^\\v\\l > 0, and by continuity there exists < ao < 1 such 
that I^^^{v°'") = 0. Therefore 

5'^,-y(y°") > d{uj). (22) 

o 1 

Consider now —S^^^{v'') = a\\d,v\\l-^\viO)\^- a^P-^^/^\\vfpXl. Since 

p ^ 5 and Q-yi^v) = 0, we have for a G [0, 1] 

^S^,.i^-) ^ aQ,{v) - |(1 - a)\vm' = -|(1 - a)|^0)p 
d 

and thus -g-S^^.y{v") ^ for all a G [0, 1], which leads to S^^^{v) ^ S^^^i^v""^). 
It follows by (22) that S^^^{v) ^ d{uj), which concludes d_^ = d{uj). 

Step 2. We construct a sequence of initial data (/?° satisfying the following 
properties : 

S.,,{^Z) < d{u),I^,,{^:) < and Q,{^Z) < 0- 
These properties are invariant under the flow of (2). 
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For a > 0, we define (/J^ by (p"{x) = a^/'^(p^^^{ax). Since p ^ 5, 7 < and 
Q"/{Viu,'y) = 0, easy computations permit to obtain 

d d 
^'S'<^,7(<^S)|"=i < 0' ^4,7(<^2)|a=i < and — Q^((^2)|a=i < 0^ 

and thus for any a > 1 close enougli to 1 we liave 

S^,,{^:) < LA^Z) < and Q.ivZ) < 0. (23) 



Now fix a a > 1 such that (23) is satisfied, and let u"'{t,x) be the solution of 
(2) with m"(0) = ifZ- Since (^2 radial, M"(t) is also radial for alH > (see 
Remark 2). We claim that the properties described in (23) are invariant under 
the flow of (2). Indeed, since from (5) and (6) we have for alH > 

S.,^{u'^{t)) = S^,^{^Z) < '5.,7(<^-,7). (24) 

we infer that I^^^{u'^{t)) 7^ for any t ^ 0, and by continuity we have 
4,7(M"(i)) < for all t ^ 0. It follows that Q^iu'^it)) ^ for any t ^ 
(if not M"(t) G ^ and thus S^^^{u'^{t)) > S^^^{ip^^-y) which contradicts (24)), 
and by continuity we have Q^{u°'(t)) < for alH ^ 0. 



Step 3. We prove that Q^{u°') stays negative and away from for all t ^ 0. 

Let t > be arbitrary chosen, define v = u°'{t) and for (3 > let be such 
that v^{x) = v{j3x). Then we have 

thus lim^_>+oo Q7(f ^) = +00, and by continuity there exists /3o such that 
Q-yiyf^^) = 0. If /(^,^(f^°) ^ 0, we keep Po unchanged; otherwise, we replace it 
by Po such that 1 < Pq < Pq, /(^,^(f^°) = and Q^i^v'^'^) ^ 0. Thus in any case 
we have 5(^,^(1'^°) ^ d{uj). Now, we have 



from the expression of and Pq > 1 it follows that 



S.,,{v) - S^,,{v^'^) ^ ^{Q,{v) - Q,{v^'^)). (25) 

Therefore, from (25), (^^(f'^") ^ and Si^^^{v^'^) ^ d{uj) we have 

Q^{v) ^ -m = 2{S^^^{v) - d{cu)) < (26) 
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where m is independent of t since S^j^-y is a conserved quantity. 



Conclusion. Finally, thanks to (26) and Proposition 6, we have 

\\xu''{t)\\l ^ -mt^ + Ct+ WxifZWl (27) 

For t large, the right member of (27) becomes negative, thus there exists 
T" < +CXD such that 

£m \\d^u''{t)\\l = +00. 
Since it is clear that — s> (f^^^j in H^(M.) when a — 1, Theorem 5 is proved. 

Remark 10 It is not hard to see that the set 

X={ve H\Ry, S^^^iv) < d{u), I^^^{v) > 0} 

is invariant under the flow of (2), and that a solution with initial data belonging 
to I is global. Thus using the minimizing character of <y9^,^ and performing an 
analysis in the same way than in [23], it is possible to find a family of initial 
data in X approaching </)^^ in (M) and such that the associated solution of 
(2) exists globally but escaped in finite time from a tubular neighborhood of 
V^w,7(see also [11,21] for an illustration of this approach on a related problem). 



6 Numerical results 



In this Section, we use numerical simulations to complement the rigorous 
theory on stability and instability of the standing waves of (2). Our approach 
here is similar to the one in [14]. In order to study stability under radial 
perturbations, we use the initial condition 

uo{x) = {l + 5p)ip^,j{x). (28) 

In order to study stability under non-radial (asymmetric) perturbations, we 
use the initial condition 

uo{x) = ipuj,^{x - 6c), (29) 

when 6c is the lateral shift of the initial condition. Since the evolution of the 
momentum for solutions of Eq. (1) is given by 



^ = -2 J \u\^VV{x)dx, (30) 

one can see that symmetry-breaking perturbations (29) do not conserve the 
momentum and thus, may give rise to drift instabilities. In some cases (when 
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the standing wave has a negative slope and the hnearized problem has two 
negative eigenvalues), we use the initial condition 

Uo{x) = (1 + Sp)ip^^y{x - (5c). (31) 

In order to demonstrate the agreement of the numerics with the rigorous sta- 
bility theory, one needs to observe that \\u — ^^1^,^11^ remains "small" in the 
case of stability but increases in the case of instability. In the latter case, how- 
ever, observing numerically that ||n — lPi^,^\\hi increases does not enable us 
to distinguish between the different types of instabilities such as total diffrac- 
tion (i.e., when limj^oo ll^lloo = 0), finite-width instability, strong instability 
or drift instability. Therefore, instead of presenting the norm, we plot the 
dynamics of the maximal amplitude of the solution and of the location of the 
maximal amplitude. Together, these two quantities give a more informative 
description of the dynamics, while also showing whether the soliton is stable. 

6.1 Stability in H}^^{M.) 

6.1.1 Strength of radial stability 

When 7 > 0, the standing waves are known to be stable in H^^^(R) for 
1 < p < 5. The rigorous theory, however, does not address the issue of the 
strength of radial stability. This issue is of most interest in the case p = 5, 
which is unstable when 7 = 0. 

For (5p > 0, it is useful to define 

r/r \ f max^ |M(x,t)| - max^.V9<^,^l 

F{6p) = max <^ \ (32 

t>o [ max^v9^,^ J 

as a measure of the strength of radial stability. Figure 2 shows the normalized 
values maXj; Iwl/max^. (f^^^^ as a function of t, for the initial condition (28) with 
uj = 4 and 7 = 1. When p = 3, a perturbation of 6p = 0.01 induces small 
oscillations and F(O.Ol) = 1.9%. Therefore, roughly speaking, a 1% pertur- 
bation of the initial condition leads to a maximal deviation of 2%. A larger 
perturbation of 6p = 0.08 causes the magnitude of the oscillations to increase 
approximately by the same ratio, so that -F(0.08) = 15%. Using the same 
perturbations with p = 5, however, leads to significantly larger deviations. 
Thus, -F(O.Ol) = 8.8%, i.e., more than 4 times bigger than for p = 3, and 
F(0.08) = 122%, i.e., more than 8 times than for p = 3. 

In [14,15], Fibich, Sivan and Weinstein observed that the strength of radial 
stability is related to the magnitude of slope (9^| |</:'a;,7l li; so that the larger 
f^Lul lv^a;,7l I2) the "more stable" the solution is. Indeed, numerically we found 
that when a; = 4, 9^^] |v5a;,7l I2 is equal to 1.0 for = 3 and 0.056 for p = 5. 
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Since when 7 = 0, the slope is positive for p < 5 but zero for p = 5, for 7 > 
the slope is smaller in the critical case than in the subcritical case. Therefore, 
we make the following informal observation: 

Observation 11 Radial stability of the standing waves of (2) with 7 > 
"weaker" in the critical case p = 5 than in the subcritical case p < 5. 

Clearly, this difference would be more dramatic at smaller (positive) values 
of 7. Indeed, if in the simulation of Figure 2 with 6p = 0.01 we reduce 7 from 
1 to 0.5 and then to 0.1, this has almost no effect when p = 3, where the value 
of F slightly increases from 1.9% to 2.1% and to 2.5%, respectively, see Fig- 
ure 3a. However, if we repeat the same simulations with p = 5, then reducing 
the value of 7 has a much larger effect, see Figure 3b, where F increases from 
8.9% for 7 = 1 to 24% for 7 = 0.5. Moreover, when we further reduced 7 
to 0.1, the solution seems to undergo collapseQ This implies that when p = 5 
and 7 > 0, the standing wave is stable, yet it can collapse under a sufficiently 
large perturbation. 



6.1.2 Characterization of radial instability for 3 < p < 5 and 7 < 

We consider the subcritical repulsive case p = 4 and 7 = —1. In this case, 
there is threshold uj2 such that (f^^^^ is stable for u > uj2 and unstable for 
uj < UJ2. By numerical calculation we found that uj2{p = 4,7 = —1) ^ 0.82. 
Accordingly, we chose two representative values of u: u = 0.5 in the unstable 
regime, and ci; = 2 in the stable regime. 

Figure 4a demonstrates the stability for = 2. Indeed, reducing the per- 
turbation from 6p = 0.005 to 0.001 results in reduction of the relative mag- 
nitude of the oscillations by roughly five times, from F(0.005) ~ 10% to 
F(O.OOl) ~ 2%. The dynamics in the unstable case u = 0.5 is also oscillatory, 
see Figure 4b. However, in this case F(0.005) = 79%, i.e., eight times larger 
than for tu = 2. More importantly, unlike the stable perturbation of 

6p = 0.001 does not result in a reduction of the relative magnitude of the os- 
cillations by ~ 5. In fact, the relative magnitude of the oscillations decreases 
only to F(O.OOl) = 66%. 

In the homogeneous NLS, unstable standing waves perturbed with Sp > 
always undergo collapse. Since, however, for p = 4 it is impossible to have 
collapse, an interesting question is the nature of the instability in the unstable 
region u < uj2. In Figure 4b we already saw that max\u{x,t) \ undergoes oscil- 
lations. In order to better understand the nature of this unstable oscillatory 

^ Clearly, one cannot use numerics to determine that a solution becomes singular, 
as it is always possible that collapse would be arrested at some higher focusing 
levels. 
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dynamics, we plot in Figure 5 the spatial profile of t)| at various values 
oft. In addition, at each t we plot where uj*{t) is determined from 

the relation 

max0^*(j) .y(x) = max|u(x,t)|. 

Since the two curves are nearly indistinguishable (especially in the central 
region), this shows that the unstable dynamics corresponds to "movement 
along the curve 0a;*(i)"- 

In Figure 6 we see that uj*{t) undergoes oscillations, in accordance with the 
oscillations of max^. \u\. Furthermore, as one may expect, collapse is arrested 
only when uj*(t) reaches a value (~ 2.86) which is in the stability region (i.e., 
above 1^2) • 

Observation 12 When 7 < and 3 < p < 5, the instability in H^^^{M.) is 
a "finite width instability" , i.e., the solution narrows down along the curve 
4'uj*{t),-y until it "reaches" a finite width in the stable region uj > 002, at which 
point collapse is arrested. 

Note that this behavior was already observed in [14], Fig 19. Therefore, more 
generally, we conjecture that 

Observation 13 When the slope is negative {i.e., di^\\(fi^^^\\2 < ), then the 
symmetric perturbation (28) with < 5p <^ 1 leads to a finite-width insta- 
bility in the subcritical case, and to a finite-time collapse in the critical and 
supercritical cases. 

6.1.3 Supercritical case {p > 5) 

We recall that when 7 > and p > 5, the standing wave is stable for 
7^/4 < u < uji and unstable for ui < uj. When 7 < and p > 5 the 
standing wave is strongly unstable under radial perturbations for any uo, i.e., 
an infinitesimal perturbation can lead to collapse. 

Figure 7 shows the behavior of perturbed solutions for p = 6 and = 1. As 
predicted by the theory, when 5p = 0.001, the solution blows up for 7 = — 1 
and 7 = 0, but undergoes small oscillations (i.e., is stable) for 7 = 1. Indeed, 
we found numerically that ci;i(p = 6, 7 = 1) ~ 2.9, so that the standing wave is 
stable for u = 1. However, when we increase the perturbation to 6p = 0.1, the 
solution with 7 = 1 also seems to undergo collapse. This implies that when 
p > 5, 7 > and u < ui the standing wave is stable, yet it can collapse under 
a sufficiently large perturbation. In order to find the type of instability for 
7 > and UJ > Ui, we solve the NLS (2) with p = 6, 7 = 1 and a; = 4. In this 
case, 6p = 0.001 seems to lead to collapse, see Figure 8, suggesting a strong 
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instability for p > 5, 7 > and uj > uji. Therefore, we make the following 
informal observation: 

Observation 14 If a standing wave of (2) with p > 5 is unstable in H^^^(R), 
then the instability is strong. 



6.2 Stability under non-radial perturbations 
6.2.1 Stability for \ < p <b and 7 > 

Figure 9 shows the evolution of the solution when p = 3,7 = l,u; = l and 
5c = 0.1. The peak of the solution moves back towards a; = very quickly 
(around t ^ 0.003) and stays there at later times. Subsequently, the so- 
lution converges to the bound state 0a;* =0.995- This convergence starts near 
X = and spreads sideways, accompanied by radiation of the excess power 
- 1 10a;* =0.995 1 12 = 2.00 - 1.99 = 0.01. In Fig 10 we repeat this simu- 
lation with a larger shift of be = 0.5. The overall dynamics is similar: The 
solution peak moves back to x = 0, and the solution converges (from the 
center outwards) to 0a;* =0.905 • In this case, it takes longer for the maximum 
to return to a; = (at t ~ 0.11), and more power is radiated in the process 
(lkoll2-|l0c.*=o.9O5|l2 = 2.00-1.81 = 0.19. We verified that the "non-smooth" 
profiles (e.g., at t = 0.2) are not numerical artifacts. 



6. 2. 2 Drift instability for 1 < p < 3 and 7 < 

Figure 11 shows the evolution of the solution for p = 3,7 = —l,u! = 1 and 
6c = 0.1. Unlike the attractive case with the same parameters (Figure 9), as 
a result of this small initial shift to the right, nearly all the power flows from 
the left side of the defect (x < 0) to the right side (x > 0), see Figure 12a, 
so that by t ~ 3, ~ 90% of the power is in the right side. Subsequently, the 
right component moves to the right at a constant speed (see Fig 12b) while 
assuming the seek profile of the homogeneous NLS bound state (see Fig 11 at 
t=8); the left component also drifts away from the defect. 

We thus see that 

Observation 15 When 1 < p < 3, the standing waves are stable under shifts 
in the attractive case, but undergo a drift instability away from the defect in 
the repulsive case. 

We note that a similar behavior was observed in the subcritical NLS with a 
periodic nonlinearity, see [14], Section 5.1. 
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6.2.3 Drift and finite-width instability for 3 < p < 5 and 7 < 

In Figure 4b, Figure 5, and Figure 6 we saw that when p = 4, 7 = — 1, 
uj = 0.5, and Sp = 0.005, the solution undergoes a finite-width instabihty 
in H^^^{W). In Figures 13 and 14 we show the dynamics (in H^{M.)) when 
we add a small shift of 6c = 0.1. In this case, the (larger) right component 
undergoes a combination of a drift instability and a finite-width instability, 
whereas the (smaller) left component undergoes a drift instability. Therefore, 
we make the following observation 

Observation 16 When 3 < p < 5, 7^/4 < u < uj2 and 7 < 0, the standing 
waves undergo a combined drift and finite-width instability. 



6.2.4 Drift and strong instability for 5 < p and 7 < 

In Figures 15 and 16 we show the solution of the NLS (2) with p = 6, 7 = — 1 
and ui = 1, for the initial condition (31) with Sc = 0.2 and Sp = 0.001. As 
predicted by the theory, this strongly unstable solution undergoes collapse. 
Note, however, that, in parallel, the solution also undergoes a drift instability. 
We thus see that 

Observation 17 In the critical and supercritical repulsive case, the standing 
waves collapse while undergoing a drift instability away from the defect. 

Note that a similar behavior was observed in [14], Section 5.2. 



6.3 Numerical Methods 

We solve the NLS (2) using fourth-order finite differences in x and second- 
order implicit Crack-Nicholson scheme in time. Clearly, the main question is 
how to discretize the delta potential at x = 0. Recall that in continuous case 

\im^dxu{x) — hm dxu{x) = —7^(0). 

Discretizing this relation with 0{h'^) accuracy gives 

u{2h) - Au{h) + 3m(0) -u{-2h) + Au{-h) - 3m(0) 



2h 2h 



-7u(0), 



when h is the spatial grid size. By rearrangement of the terms we get the 
equation 

- u{2h) + Au{h) + [2/17 - 6]m(0) + Au{-h) - u{-2h) = 0. (33) 
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When we simulate symmetric perturbations (section 6.1), we enforce symme- 
try by solving only on half space [0, +00). In this case, because of the symmetry 
condition u{—x) = u{x), (33) becomes 

[2/17 - 6]n(0) + 8u{h) - 2u{2h) = 0. 



A Proof of Lemma 10 



The proof for Ll^ being similar to the one of Lj^ we only deal with 
Lj^^. The form Bj^^ can be decomposed into Bj,^ = BJi + -B7,2,w with 
Bl^ : H\R) X H^R) M and S7^2,a. : L'^W x ^^(K) ^ ^ defined by 

Bi,2,c,{'f, z) = w /k vzdx - j^p^l~^vzdx. 

If we denote by Ti (resp. T2) the self-adjoint operator on L^(]R) associated 
with 57,1 (resp 57,2,^.), it is clear that D{T2) = /.^(R) and 

D{Ll^)=D{T,). 

If we take v e H'^(R) such that v{0) = 0, and put w = -dlv e ^^(M), it 
follows that for any z G H^{R) we have 

B^i{v,z) = / dxvdxzdx = {w,z)2. 

Thus V G D{Ti), and we can deduce that Ti is a self-adjoint extension of the 
operator T defined by 

T = -dl D{T) = {v G H'iR); t;(0) = 0}. 

On the other hand, using the theory of self-adjoint extensions of symmetric 
operators, one can see (see [2, Theorem 1-3.1.1]) that there exists a G M such 
that 

D(Ti) = {v G H\R) n H\R \ {0}); dX^^) - dMO~) = -av{0)}. 
Now, take v G D{Ti) with v{0) ^ 0. Then 
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/O r+oc 
{-dlv)vdx+ / {-dlv)vdx 
-oo Jo 



/O r+oc 
\d^v\'^dx + v{0)d^v{0+) + / \d^v\^dx 
-oo Jo 

\d^v\^dx - av{Of 

which should be equal to 

Bj,iv,v) = [ \d,v\^dx-^v{Of. 
Jr 

Thus 7 = a, and the lemma is proved. 



B Proofs of Lemmas 11 and 13 



Proof of Lemma 11. We start by showing that (i) and (iii) are satisfied. We 
work on Ll^^ and L^i^. The essential spectrum of Ti (see the proof of Lemma 
10) is acssiTi) = [0, +oo). This is standard when 7 = and a proof for 7 7^ 
can be found in [2, Theorem 1-3.1.4]. From Weyl's theorem (see [30, Theorem 
lV-5.35]), the essential spectrum of both operators Lj^^ and Ll^^ is [u, +00). 
Since both operators are bounded from below, there can be only finitely many 
isolated eigenvalues (of finite multiplicity) in {—oo,uj') for any uj' < uo. Then 
(i) and (iii) follow easily. 

Next, we consider (ii). Since (p^^^^ satisfies Ll^i^ip^^^ = and ip^^^ > 0, the first 
eigenvalue of L2 is and the rest of the spectrum is positive . This is classical 
for 7 = and can be easily proved for 7 7^ 0, see [5, Chapter 2, Section 2.3, 
Paragraph 3]. Thus to ensure that the kernel of is reduced to spa.n{iip^^^} 
it is enough to prove that the kernel of Lj^ is {0}. It is equivalent to prove 
that is the unique solution of 

Ll^u = 0,ueD{LjJ. (B.l) 

To be more precise, the solutions of (B.l) satisfy 



u e H\R\{0})nH\R), (B.2) 



-dlu + LJU- p^l~^u = 0, (B.3) 
9,.m(0+) - 9,m(0-) = -7m(0). (B.l) 

Consider first (B.3) on (0, +00). If we look at (3) only on (0, +00), we see that 
ip^^^ satisfies 

- 5^V5^,7 + uj^Puj,'y - V^S,^ = on (0, +00). (B.5) 
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If we differentiate (B.5) with respect to x (which is possible because ip^ia 
is smooth on (0, +00)), we see that d^^^^^^ satisfies (B.3) on (0,+cx)). Since 
we look for solutions in L^(M) (in fact solutions going to at infinity), it is 
standard that every solution of (B.3) in (0, +cx)) is of the form [idx'^^^^, /i G M 
(see, for example, [5, Chapter 2, Theorem 3.3]). A similar argument can be 
applied to (B.3) on (—00, 0), thus every solution of (B.3) in (—00, 0) is of the 
form vdx^uj^-yi G IR- 

Now, let u be a solution of (B.2)-(B.4). Then there exists /x G M and z/ G M 
such that 



u = vdx'^uja on (-00, 0), 
u = fidx(Puj,'y on (0, +cxd). 

Since u G H^(R), u is continuous at 0, thus we must have /i = — z/, that is u 
is of the form 



u = -iJ.dx(puj,-y on (-CX),0), 

U = l2da:(Puj,-y OU (0, +Oo), 



-/i 



m(0) = -IJ,dx(Puj,-y{0-) = l^dxiPuj,-y{0 + ) = -^lVuj,'y{0). 



Furthermore, u should satisfies the jump condition (B.4). Since ip^^^ satisfies 
if we suppose /i 7^ then (B.4) reduces to 

But from (4) we know that 

It is a contradiction, therefore /i = 0. In conclusion, m = on M, and the 
lemma is proved. 



Proof of Lemma 13. We recall that VI is defined with the help of a bilinear 
form Bl (see (9)). To prove the holomorphicity of (Lj) it is enough to prove 
that (-Bl) is bounded from below and closed, and that for any v G H^{M) the 
function Bj{v) : 7 ^ Bj{v,v) is holomorphic (see [30, Theorem VII-4.2]). It is 
clear that is bounded from below and closed on the same domain H^(M.) for 
all 7, thus we just have to check the holomorphicity of Bl{y) : 7 t-^ -Bi (f, v) 
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for any v G if^(M). We recall the decomposition of into i?7i and BI2 (see 
(A.l)). We see that ^7 i{v) is clearly holomorphic in 7. From the explicit form 
of ip^ (see (4)) it is clear that 7 1— > ^p^~^{x) is holomorphic in 7 for any x G M. 
It then also follows that 7 ^ Bi 2{v) is holomorphic. □ 

Remark 18 There exists another way to show that (Lj) is a real-holomorphic 
family with respect to 7 G M. We can use the explicit resolvent formula in [2], 

(Ti - = i-dl - kY' + 27M-n + 2A;)-1(G;(0, ■)Gk{-), 

where k"^ G p(T'i), Imfc > 0, Gk{x) = {i/2k)e'^^^^\ to verify the holomorphicity. 

C Proof of Lemma 20 

First, we indicate how the extension of / to (— cxd,+oo) can be done. We 
see by the proof in [38, Theorem XII. 8] that the functions 7(7) and A(7) 
defined in Lemma 15 exist, are holomorphic and represent an eigenvector and 
an eigenvalue for all 7 G M, since {Lj) is a real-holomorphic family in 7 G M. 
Namely we can repeat the argument of Lemma 15 at each point 7 and on each 
neighborhood of 7. This is possible because the set {(7, A); 7 G M, AG p(L]')} 
is open and the function (L]^— A)~^ defined on this set is a holomorphic function 
of two variables ([38, Theorem XII. 7]). 

Secondly, as it was observed in [12,14], the eigenvectors of VI are even or odd. 
Indeed, let ^ be an eigenvalue of L\ with eigenvector v G D{LJ). Then clearly 
V with v{x) = v{—x) is also an eigenvector associated to ^. In particular, v 
and V satisfy both 

-d^v + {00 - - pip'y^v = on [0, +00), 

thus there exists 77 G M such that v = rjv on [0, +00) (this is standard, see, for 
example, [5, Chapter 2, Theorem 3.3]). If u(0) 7^ 0, it is immediate that 77 = 1. 
If w(0) = 0, then dxv{Q+) 7^ (otherwise the Cauchy-Lipschitz Theorem leads 
to f = 0), and it is also immediate that 77 = —1. Arguing in a same way on 
(—00, 0], we conclude that v is even or odd, and in particular v is even if and 
only if v{Q) ^ 0. 

Finally, we prove the last statement only for the case 7 < since 
the case 7 > is similar. We remark that d^^o is odd. Since 
lim^^o(/(7),^xV5o)2 = ||5x<Po||2 7^ 0' have (/(7), <9x</3o)2 7^ for 7 close 
to 0, thus /(7) cannot be even, and therefore /(7) is odd. Let 700 be 

7oo = inf{7 < 0; 7(7) is odd for any 7 G (7, 0]}. 
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We suppose that 700 > —00. If /(700) is odd, by continuity in 7 of /(7) with 
value, there exists e > such that fi'joo — ^) is odd which is a contradiction 
with the definition of 700, thus /(700) is even. Now, /(700) is the hmit of odd 
functions, thus is odd. The only possibility to have /(700) both even and odd 
is /(7oo) = 0, which is impossible because /(700) is an eigenvector. □ 



D Proof of Proposition 6 



For a G N \ {0}, we define V°'{x) = ■yae It is clear that JrV°'{x) = 7 

and V"" 75 weak-Tir in H^^{M.) when a — > +00. 

We begin by the construction of approximate solutions : for 

' ' ' (D.l) 

m(0) = Mo, 

and thanks to [7, Proposition 6.4.1], for every a G N \ {0} there exists > 
and a unique maximal solution u" G C([0, T'^), i/^(M)) n CH[0, T"), if-i(M)) 
of (D.l) which satisfies for all t G [Q.T") 



E%u%t))=E^{uo), (D.2) 
\K{t)h = \\uo\\2, (D.3) 

1 If 1 

where E°'{v) = -\\dxv\\l — / Prfx Moreover, the function 

2 2 Jr p + 1 ^ 

: t ^ J^x^\u''{t,x)\'^dx is by [7, Proposition 6.4.2], and 



dtr = ^lm J U^xd^u'^dx, (D.4) 

d^,r = SQt^{u'^) (D.5) 

where is defined for v G if^(]R) by 

= \\dM\l + \ lx{d^Vn\v\'dx - ^^^\\v\C\. 



Then, we find estimates on Let M ^ ||^io||//i(R) (an exact value of M will 
be precise later). We define 

t" = sup{t > 0; ||u"(s)||hi(r) ^ 2M for all s G [0,t)}. (D.6) 
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Since u"" satisfies (D.l), we have 

sup ||9tM"||Loc([0,ta)^^-l(IR)) ^ C, 

aGN\{0} 

and thus for all t G [0, and for all a G N \ {0} we get 

h^it) -uo\\l = 2 f\u%s) - no, dtu%s))2ds ^ Ct (D.7) 
Jo 

where C depends only on M. Now we have 
1 

p + 



Ti\hX+i-\ho\[pXi)= / {u''-uo)\su'' + ii-s)uo\pdxds 



which combined with Holder inequality, Sobolev embeddings, (D.6) and (D.7) 
gives 

^(IkllSl - hoRl) ^ Ct'l\ (D.8) 

Moreover, using (D.6), Sobolev embeddings, Gaghardo-Nirenberg inequality 
and (D.7) we obtain 



v^'^(Kp-KP; 



^ Ct^l\ (D.9) 



Combining (D.2), (D.3), (D.8) and (D.9) leads to 
l|w''(^)l|j/i(R) ^ + Cit^l^ + t^/^) for all t G [0,^) and for all a G N \ {0}, 

and choosing Tm (depending only on M) such that C(tI(^ + t2/^) = 3M^ we 
obtain for all a G N \ {0} the estimates 

||^''IU°°([o,Tm),hi(m)) ^ 2M, 

\\^tU''\\Loo{[O^TM),H-^R)) ^ C. 

In particular, it follows from (D.IO) that Tm < for all a G N \ {0}. 

Now we can pass to the limit : thanks to (D.IO) there exists 
u G L~([0,TM),i7i(M)) such that for all t G [0,Tm) we have 

u^it) u{t) weakly in H\R) when a +oo, (D.ll) 

which immediately induces that when a +oo, 

\u''{t)\P^\^{t) \u{t)\P~\{t) weakly in H~\R). (D.12) 

In particular, thanks to Sobolev embeddings, we have 

u°'{t,x) u{t,x) a.e. and uniformly on the compact sets of M, 



33 



and it is not hard to see that it permit to show 

W ^ M7(5 weak-* in (D.13) 

Since satisfies (D.l), it follows from (D.ll), (D.12) and (D.13) that u sat- 
isfies (2). Finally, by (6) and (D.3), we have 

n'^-.uinC([0,TM),L2(M)), 

thus, from Gagliardo-Nirenberg inequality and (D.IO), we have 

u^^umC{[0,TM),LP+\R)), 

and by (5) and (D.2) it follows that 

u"" ^uin C([0, Tm), H\R)). (D.14) 

We have to prove that the time interval [0,Tm) can be extended as large as 
we need. Let < T < T^^ and 

M = sup{||M(t)||Hi(E),te [0,T]}. 

If Tm ^ T, there is nothing left to do, thus we suppose Tm < T. From (D.14) 
we have ||tt"(7M)||_ffi(R) ^ M for a large enough. By performing a shift of time 
of length Tm in (2) and (D.l) and repeating the first steps of the proof we 
obtain 

u''^umC{[TM,2TM),H\R)). 

Now we reiterate this procedure a finite number of times until we covered the 
interval [0,T] to obtain 

^ M in C([0,T],i/i(M)). (D.15) 

To conclude, we remark that (7) follows from the same proof than [7, Lemma 
6.4.3] (computing with ||e^'^'' a;u(t)||2 and passing to the limit e 0), thus we 
have 

11^^(^)112 = ll^^o||2 + 4 / Im / u{s)xdxu{s)dxds. (D.16) 
From (D.4), Cauchy-Schwartz inequality and (D.IO) we have 

dt{\\xu''mt) ^C\\xu'^{t)U 

which implies that 
Since in addition we have 

xu"'{t,x) -^xu{t,x) a.e.. 
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we infer that 

xu°'{t,x) xu{t,x) weakly in L^(]R). 

Recalling that 

dxu"' — > dxU strongly in L^(]R) 
we can pass to the limit in (D.16) to have 

\\xu^{t)\\2 \\xu{t)\\2. 

On the other side, since we have (D.5) and (D.15), we get (8). 

Remark 19 Our method of approximation is inspired of the one developed in 
[9] by Cazenave and Weissler to prove the local well-posedness of the Cauchy 
problem for nonlinear Schrodinger equations. Actually, slight modifications 
in our proof of Proposition 6 would permit to give an alternative proof of 
Proposition 2. 
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Fig. 2. maxx\u\/maxx function of t for = 4, 7 

line) and 5p = 0.08 (solid line), (a) p = 3 (b) p = 5. 
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Fig. 3. maxx\u\/maxx'^uj^'y as a function of t for w = 4, 5^ = 0.01, and 7 = 1 (solid 
line), 7 = 0.5 (dashed line) and 7 = 0.1 (dots), (a) p = 3 (b) p = 5. 
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Fig. 4. maxa;|u|/ma3;2:(/3a;,7 as a function of t for p = 4, 7 = —1, 5p = 0.001 (dashed 
line) and 6p = 0.005 (solid line), (a) w = 2; (b) a; = 0.5. 
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Fig. 5. |n(x,t)| (solid line) and (puj^:(t){x) (dots) as a function of x for the simula- 
tion of Fig. 4b with 6p = 0.005. (a) t = {oj* = 0.508) (b) t = 9 {oj* = 1.27) 
(c) t = 10.69 {oj* = 2.86) (d) t = 12 {u* = 1.43) (e) t = 15 (w* = 0.706) (f) 
t = 20 (lo* = 0.58). 
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Fig. 6. LO* as a function of t for the simulation of Fig 5. 
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Fig. 7. maxa; \u{x,t)\/ maxx (p^^-y as a function of t for p = 6,uj = 1 and 7 = — 1 
(dashed line), 7 = (dots), 7 = +1 (solid line), (a) 6p = 0.001 (b) 6p = 0.1. 
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Fig. 8. maxx\u{x,t)\/ maxx 'fuj.'y as a function of t for p = 6,u; = 4, 7 = 1 and 
6p = 0.001. 




Fig. 9. (solid line) and (puj*=o.995ix) (dashed line) as a function of x. Here, 

p = 3,u; = l,7 = l and 6c = 0.1. 
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Fig. 12. (a) The normalized powers \u\'^dx/ \uo\'^dx (solid line) and 

I-oo \'^\'^ / l^oprfa; (dashed line), and (b) location of m.sxo<,x\u{x,t)\ (solid 
line) and of maXa;<o |u(x,t)| (dashed line), for the simulation of Figure 11. 
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Fig. 13. u{x,t) as a function of x. Here p = 4, 7 = —1, co = 0.5, 6p = 0.005, and 
6c = 0.1 
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Fig. 14. (a) The value, and (b) the location, of the right peak maxo<a; (solid 
line) and left peak max^<o i)| (dashed line), for the simulation of Figure 13. 
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Fig. 15. |n(x,t)| as a function of x, at various values of t. Here, p = 6, 7 = — 1, 
uj = l,5c = 0.2 and 5p = 0.001. 
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Fig. 16. (a) max^; 1^(2;, max^ (^^^^ (b) location of max^; |n(x, t)| and (c) The 
normalized powers \u\'^dx/ |nopdx (solid line) and \u\'^dx/ |tiop(ix 
(dashed line), for the solution of Fig. 15. 
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